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Abstract 

We construct a classical dynamical system whose phase space is a certain infinite- 
dimensional Grassmannian manifold, and propose that it is equivalent to the large N 
limit of two-dimensional QCD with an 0{2N + 1) gauge group. In this theory, we find 
that baryon number is a topological quantity that is conserved only modulo 2. We also 
relate this theory to the master field approach to matrix models. 
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1 Introduction 



An interesting type of a theory of the strong interactions is one in which there is a 
topological quantity which can be related to hadrons; several theories relate such a 
quantity to baryon number. The plausibility of such theories may be enhanced by 
their being consistent with the well-accepted gauge theory of the strong interactions 
- i.e. with QCD. 

One paper by Rajeev proposes such a theory in two dimensions, and shows its 
equivalence to the large N limit of SU(N) QCD pQ. The theory is a classical 
one in which the phase space is an infinite dimensional Grassmannian manifold, 
and baryon number is the topological invariant that corresponds to the Z-many 
connected components of the manifold. Rajeev's theory is tied into the gauge theory 
of the strong interactions via the equation for the mass of the mesons, which agrees 
with the one derived by 't Hooft in the large N limit of SU(N) gauge theory [2]. 

Here, we propose an analog of Rajeev's theory for the large N limit of 0(2N+ 1) 
QCD in two dimensions^, and arrive at an interesting consequence regarding baryon 
number non-conservation. For our phase space, we construct a different infinite 
dimensional Grassmannian manifold for which the meson equation is the same as 
the one arrived at in the 0(2N + 1) gauge theory. Our phase space has a topological 
invariant of its own, which also corresponds to baryon number. Unlike the SU(N) 
case, in this case there are only two connected components; particles with even 
baryon number are assigned to one component, and those with odd baryon number 
to the other component. We thereby discover that baryon number is conserved only 
modulo two, and baryons may annihilate or be created in pairs. This fact, as we 
will explain, is independently true of QCD with an 0{2N + 1) gauge group. 
1 A different analog was considered in [3l H] . 

1 



We also show that our theory is related to the master field approach to large N 
matrix models via a master field whose commutation relations match those of our 
theory. 

This paper is organized as follows. In Section [2] we review 't Hooft's derivation 
of the equation for mesons in two dimensions in the context of the large N limit of 
U(N) QCD. In Section [3], we review the theory developed by Rajeev, concentrating 
on the derivation of the meson equation and the topological invariant corresponding 
to baryons. In Section H] we present our analog for QCD with an orthogonal gauge 
group along with the topological invariant corresponding to baryons. In Section [5] 
we show that our result about baryon number non-conservation agrees with two- 
dimensional QCD with an odd orthogonal gauge group. In Section [6] we discuss 
the relation between our model and the master field approach to matrix models for 
large N QCD. We conclude in Section [7] with several suggestions for further work. 

2 Planar diagrams and mesons in two dimensions 

The planar diagram theory developed by 't Hooft [5] and elaborated upon by Witten 
|6j provided a simplification of QCD gauge theory, 't Hooft considered QCD with 
color gauge group SU(N) in the limit N — > oo with g 2 N held fixed, and arrived at 
the notion that in this limit, only planar diagrams need to be considered, all others 
being suppressed by factors of 1/N. Still, however, calculations remain complicated 
even when they include only the planar diagrams. 

For the case of two-dimensional QCD, 't Hooft showed that a further simplifica- 
tion arises, and derived an equation for the meson spectrum [2J. The derivation is 
briefly reviewed as follows: 
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We start with the QCD lagrangian 

C = - ^(i^D^ + m)^ (1) 

where F^ v is the field strength, the covariant derivative, ip the quark wave 
function, and m the quark mass. In two dimensions, the index [i runs over {0, 1}. 
Switching to light-cone coordinates, where x ± = (x 1 ±x°)/\^2, p± = (pi ±p )/\/2, 
A± = (Ai ± A ) I y/2, and g a b — — 1, we impose the gauge A_ = A + = and the 
lagrangian simplifies to: 

£ = -^Tr(<9_A+) 2 - ^(zyfy + m + gj-A + )^. (2) 

There is only one vertex, the gauge fields do not interact with themselves, and the 
feynman rules simplify considerably. 

To derive the equation for the meson, we first need the dressed propagator, 
G(k,m), for the quark. Let iT(k) be the amplitude for the (planar) irreducible 
quark self-energy diagram. Then G(k, m) is a sum of diagrams created from the 
irreducible self-energy blob, and is equal to 

G{K m) = m 2 + 2k+k„-k_Y{k)-ie (3) 
An expression for T(k) can be derived from a bootstrap equation from which a 
logarithmic UV divergence is removed and into which an IR cutoff A is introduced. 
The expression is 

r W = r(M = -^f^-^, 

7T \ A p-j 

and using it in equation (j3J) gives 

— 



G(k, m) 



m 2 + 2k + k^ — g 2 /n + g 2 \k^\/ix\ — ie 
—ik_ 

M 2 + 2k + k_ + g 2 \k_\/7r\ - it 
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(4) 



where M 2 = m? — g 2 /ir. Since A is small, we see that the poles of G(k) occur at 
k + —>■ oo, which means there is no physical single quark state - no free quarks. 

To find the spectrum for mesons, we consider a blob out of which come a quark 
and an antiquark. 

One can derive a bootstrap equation for such a blob. Let ip(p,r) represent a 
blob out of which go a quark with mass mi and momentum p and an antiquark of 
mass m2 and momentum p — r. Then a bootstrap equation for it is given by 

f d 2 k 1 

ip(p,r) = -j 2 A g 2 G(p - r, m 2 )G(p, rnx)^-ip(p + k, r). (5) 

Defining 0(p_,r) = J dp + if){p + ,p-,r), and substituting the expressions for G(k,m) 
from equation (J4]), we get after some algebra 

where /i 2 7r/g 2 = r ■ r = — 2r + r_ (i.e. \x is the meson mass in units of g/y/W), 
ai = irm 2 /g 2 — 1, x = p_/r_, "P" stands for the principal value integral, i.e. 

p f <t>{y) d = 1 /" + ie) rf + 1 /■ 0(y - ze) 
J (y — x) 2 2 J (y + ie — a;) 2 2 7 (y — ie — x) 2 ' 

which is finite, and where the infrared cutoff A has disappeared. Equation ([2]) is 
known as 't Hooft's integral equation for mesons. The spectrum (i.e. the eigenvalues 
yU 2 ) can be shown to be a discrete set of positive eigenvalues if mf > 0. 

When we replace U(N) with O(M) as the gauge group (where M may be even 
or odd), the derivation of the integral equation is the same so the same equation 
holds for mesons in two-dimensional QCD with an O(M) gauge group. 
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3 SU(N) QCD on a Grassmannian Manifold 

Now we review the main ingredients of Rajeev's construction of the large N limit of 
two-dimensional QCD as a classical dynamical system [I]. 

There are three necessary ingredients for a classical dynamical system: a phase 
space, a symplectic form, and a Hamiltonian. Here, the phase space is given by 
an infinite dimensional Grassmanian manifold, described in two equivalent ways. 
The first description is the following set of operators on a polarized Hilbert space 
H = H+®H-: 

Gn = {$ | $t = $ ; $ 2 = 1 ; [ e , $] is Hilbert-Schmidt}, (7) 

where e = ±1 on 7i±, and an operator A is Hilbert-Schmidt if tr(A T A) < oo. The 
Hilbert space is taken to be the space of square integrable complex valued functions 
on the real line (or the circle, with radius taken to infinity); it is spanned by Fourier 
modes: F{0) = E-oo F me ime ; H+ is defined as the span of {e imd , m > 0}, and H- 
is its orthogonal complement. 

It will be convenient at times to rewrite this description in terms of M = $ — e: 

Gr l = {M | M f = M ; M 2 + eM + Me = ; [e, M) is Hilbert-Schmidt}. (8) 

The second description of the space, equation (jSJ), comes from the following realiza- 
tion: each $ G Gr\ can be diagonalized to e via the action of a unitary transforma- 
tion g on the same Hilbert space: gcj)g^ = e, where by "unitary" we mean that g is 
an element in U res (TC) defined by 

U res (H) = {g | g ] g = 1; [e, g] is Hilbert-Schmidt} . (9) 
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The stabilizer of e under this action of the unitary group is U(TC+) x U{7i ). 
Therefore, we may view Gr\ as the coset space 

Gn = U(H + ) x U{H.) ' (10) 
which can be recognized as an infinite dimensional analog of finite dimensional 
Grassmannian spaces such as U(n)/[U(r) x U(n — r)\. 

To derive equations of motion, we need a symplectic form - which would give the 
Poisson Brackets - and a Hamiltonian. The symplectic form on this space is given 
by 

iv(U,V) = - l -Tr<S>[U($),V($)}, (11) 

where U, V are tangent vector fields to Gr\ (they, too, are hermitian operators, and 
they satisfy V($)$ + = 0). 

While we will not be using this symplectic form for our case, we review here some 
of its properties: in addition to u(U, V) being a closed, non-degenerate two-form, it 
is invariant under the unitary group action, $ — > gQg^ orM-> gMg< + g[e, g*]. For 
the infinitesimal form of this action, $^$ + l / u = $-M[-u,e + M\ where e tu = g 
and u is Hermitian, there is a function f u associated with each V u defined by 

u(V w •) = df u - 

namely 

f u {M) = - l -TrMu. 

We note here that while in a finite dimensional space, the associated function 
f u of any vector field always exists (the symplectic form, which is nondegenerate, is 
"invertible"), this is not always the case for an infinite dimensional space, but such 
a function does exist here. 
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The Poisson Brackets of two such functions are defined in terms of u by 

{/«, fv} = w(y u , v v ) . 

We get the following relation: 

{f u ,f v }(M) = f i[uM (M) - l -Tr[e,u]v . (12) 

We will now translate this into integral kernel language: the integral kernel M(x, y) 
of the operator M, (also known as a master field - see Section [6]), for any F G 7i, is 
defined by 

(MF){x) = J M(x,y)F(y)dy . (13) 

The kernel e(x, y) = e(x — y) for the operator e is defined similarly (and is known as 
the Hilbert transform operator). Furthermore, a trace such as TrMu is rewritten 
as / dx dy M (x , y)u(y , x) . Now we can rewrite equation ( fl2i) in terms of integral 
kernels: 

?-{M(x,y),M(z,w)} =8(x-w)[M + e](z,y) — 8(y — z)[M + e] (x, w) (14) 

(one can go back from this form to equation (TT2|) by multiplying by u(y, x)v(w, z) 
and integrating over all four variables.) The Hamiltonian is taken to be 

H(M) = J dxdyh(x - y)M(x, y) - ^g 2 J dxdyG(x - y)M(x, y)M(y, x) , (15) 

where the first term is the kinetic term, the second is the potential (interaction) 
term, and g is a constant parameter. The kernels h(x — y) and G(x — y) are given 
by 

% % 1 

h(x - y) = ^{-8'{x -y) + -sgn(x - y)) , G{x -y) = --\x-y\ , 



which are the kernels of the Fourier transforms of 



The kinetic term above is derived using 2-dimensional space with quasi-light cone 
coordinates: let u = x° — x 1 , x = x 1 , with metric ds 2 = du(du + 2dx). The momenta 
in these coordinates are p u = p (associated with energy) and p — p x — Po + Pi- In 
these coordinates, the invariant mass ji 2 is given by 



= / dz [h(x - z)M(z, y) - M(x, z)h(z - y)] 

+g 2 J dz G(y - z)e(x, z)M(z, y) - G(z - x)e(z, y)M(x, z) 
+g 2 J dz M(x, z)M(z, y)[G{y - z) - G(z - x)] . 

To arrive at the meson equation, we take a linear approximation around the vacuum 
$ = e (M = 0), so we neglect the terms quadratic in M: 



2pp u - p 2 = n 2 , 



so 



1 . a 2 . 



where p u can be interpreted as the kinetic energy. 

Now, the equation of motion, with t as the time variable, is 



1 dM(x,y,t) 

2 Ft 



= l -{H,M(x,y)} 




= dz[h(x - z)M(z,y) - M(x,z)h(z - y)\ 



+g 2 / dzG(y - z)e(x,z)M(z,y) - G(z - x)e(z,y)M(x,z) . (16) 
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Translating to momentum space, where 

M(p,q) = J dxdye i{ - px+(iy) M{x,y) , (17) 

the linearized equation becomes 

i dM(v at) f dv 

-[h(p)-h(q)]M(p,q) + g 2 (sgn(p) - sgn(q)) \ —M(p-r,q-r) . (18) 



2 dt 1 ^ ' y " y ' v x ' v "J r 

Now, define 

P = p-q, £ = p/P, x(P,€) = PM{P(,-(l-£)P) . (19) 

The quantity P has the meaning of momentum, since the translation M(x, y) — > 
M(x+a, y+a) yields in momentum representation M(p, q) — > e^~ pa+ga ^ M(p, q). Fur- 
thermore, the constraint (see equation (jSJ)) M 2 + Me + eM = 0, which to first order 
in M is Me + eM = or in kernel language / dy M(x, y)e(y, z) + e(x, y)M(y, z) = 0, 
has the Fourier transform 

(sgn(p) + sgn(q))M(p, q) = . 

We also have M^ = M, or M*(y,x) = M(x,y), which translates to M*(p,q) = 
M(q,p) in momentum space. Putting these together we note that M(p,q) may be 
non-zero only when the signs of p and q are opposite, and due to the hermiticity 
condition, it is sufficient to consider the case where p > 0, q < 0. These conditions 
also imply that £ of equation (fl9l) ranges between and 1. Now, changing from 
p, q, M(p, q) to P, £, x(P, £), and using the ansatz M(p, q, t) = e~ iPuU M(p, q) we get 

(2p u P - P 2 ) X (0 = M 2 X (0 = + y^)x(0 + V 7^pX(0 , (20) 

which is 't Hooft's integral equation for mesons (see equation (jSJ)). 

Now we turn to baryons, which are topological solitons in this classical theory, 
and discuss the relevant topological properties of Gr\. for an element $ G Gr±, 
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let V$~ be the subspace of Ti which is the +1 eigenspace of $, and let V$ be the 
subspace of 7i which is the —1 eigenspace of $ (the condition $ 2 = 1 leads to the 
fact that <3? has eigenvalues A = ±1); also, let Mq> = $ — e as before. The fact that 
[e, $] is Hilbert-Schmidt implies (see [7]) that the number of independent vectors 
in 7i which have eigenvalue +1 under the action of e but become —1 eigenvectors 
under $ is finite, i.e. dim{y§ fl 7i + ) is finite; similarly, dim(V$ fl TL-) is finite. 
Define 



The quantity J($) can take on any integer value and is known as the virtual rank, 
or index, of the operator It is a topological invariant: smoothly varying $ leaves 
the index fixed. This index divides Gr\ into Z connected components. 

The integer J($) corresponding to each $ is the topological invariant which 
corresponds to the baryon number. 



4 0(2N + 1) QCD on a(nother) Grassmannian 



Here we construct a manifold S, a real analog of Gr\, as the phase space for a 
classical dynamical system, which we propose to be equivalent to the large N limit 
of 0(2N + 1) QCD. The hamiltonian will be taken to be the analog of the one in 
Section [3l as for the symplectic form, we take it to be the canonical one on the loop 
group L02n, which is closely related to <S. We will show that mesons in this theory 



J($) = / + ($) - /_($) = dim(V$ n H+) - dim(V£ n H-) , 



(21) 



which can also be written 




(22) 
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satisfy 't Hooft's integral equation. In the next section we will discuss baryons in 
this theory. 

As in Section [31 the phase space can be described in two ways. The definition 
analogous to equation ([7]) is 

S = {$ | $ T = $ ; $ 2 = 1 ; [e, $] is Hilbert-Schmidt} , (23) 

where here $ are operators with real matrix elements, and we use the transpose 
instead of hermitian conjugate. Again, there is a description in terms of M = $ — e: 

S = {M | M T = M ; M 1 ■ (M ■ Mt 0: [e, M] is Hilbert-Schmidt} . (24) 

Following a reasoning similar to that in Section El we get the analog of equation 
(HO}: 

Ores(H) (25) 



0(H+) x 0(H-) 
where 

Ores(H) = {g\g T g = l;[e,g] is Hilbert-Schmidt} . (26) 

Turning to the symplectic form, we note that for S, the form analogous to the one 
defined in equation ffTTl) is not nondegeneratej Furthermore, the Poisson Brackets 
of equation (fl4|) are inconsistent with the constraint $ T = $, which translates into 
M(x,y) = M(y,x). Therefore, we shall describe a different space, the loop group 
L02ni which is closely related to S and on which there is a canonical symplectic 
form. 



2 Proving nondegeneracy for the form ui in equation pip (see [T]) involves looking at lo at the chosen point 
$ = e, where we have u> € {U, V) oc tr(v)v — v^u) for appropriate u, v. At this chosen point, uo c is non-degenerate, 
but for S, we would have tr(u T v — v T u) which is identically zero. Unlike here, in [4], a complexified version of 
equation P5)) was used, and their analog of equation pd|) is still non-degenerate. 
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The group L0 2n , which is the loop group of the orthogonal group 0(2n), is 
defined as the group of smooth maps from the circle to a closed curve in 0(2n): 



L0 2n = {7 : S 1 — 0{2n)} , (27) 

where the group multiplication is given by pointwise multiplication in 0(2n). Since 
the fundamental group of 0(2n) is Z 2 : 

7r 1 (0(2n)) = Z 2 , (28) 

we see by definition that L0 2n has two connected components. We shall use them 
below in connection with baryon number. 

The group L0 2n is related to S because it is embedded in O res (TC): let Tt be the 
Hilbert space of square integrable functions from the circle to R 2n o Then 7 6 L0 2n 
acts on F e Ji by pointwise multiplication given by an operator M 7 : 

(M,.F)(6) = 1 (6)F(6) . (29) 

If we write 7(0) = Elke ike , 7 fc G 0(2n), and F(6) = E^e ae , F e R 2n , then in 
terms of the standard basis for 7i, M 7 is a Z x Z matrix, with each entry M pq being 
a 2n x 2n matrix, given by M pq = j p -q. It can be shown ([7]) that this action is in 
O res (Tl), so that L0 2n is a subgroup of O res {7i). 

Another fact relating L0 2n and S is their number of connected components: it 
is a property of S that 

dim(y§ n W+) = dim(VJ|" n HJ) , (30) 



3 Note that this Hilbert space, L 2 (S X ; R 2n ), may be identified with the Hilbert space £ 2 (R; C) defined as Ti. 
at the beginning of Section 3. See [7]. 
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i.e. !+($) = /-($) = (the definitions of V$, I±($), and J($) are analogous to 
the ones corresponding to U res given in Section [3]). From this follows that J($) = 
always; however, another property of S is that smoothly varying $ may change 
/($) only by multiples of 2. This means that the parity of the dimension is 
a topological invariant; it divides S into two connected components, namely 
even and odd. This can be written as follows: 

S($) = mod 2 = -Ifr(M)| w+ mod 2 , (31) 

so H($) is either or 1, and it tells us explicitly to which connected component 
<3> G S belongs. So L0 2 n and S both have two connected components. We shall use 
this in section 

Now we construct the symplectic form u on LOm i n terms of the Killing form 
of o(2n), the Lie algebra of 0(2n). We begin by defining the Lie algebra Lo 2 n 
associated with L0 2n . It is the set of smooth maps from the circle to o(2n): 

Lo 2n = {r ] :S 1 ^ o(2n)} , (32) 

where the commutator is defined by pointwise commutators in o(2n). The Killing 
form on o(2n) - which is a symmetric, nondegenerate, invariant bilinear form - is 
defined by 

K(X, Y) = Tr(adX adY) , (33) 

where X, Y G o(2n), Tr denotes the trace of a matrix, and adX denotes the operator 
in the adjoint representation of o(2n) corresponding to X. The calculation of the 
Killing form on o(2n) is straightforward and gives 

K(X, Y) = (2n - 2)Tr(XY) . (34) 
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Now, a symplectic form on L0 2n is a two-form, namely a map which takes two 
vector fields and gives a real number. Since a Lie algebra in fact consists of vector 
fields at the identity of the group manifold, we can define a symplectic form on the 
identity component Q0 2n of the loop group as follows (see [7]): 

1 r 27r 

U ( V ,{) = — K(r,(e),e(e))d9 . (35) 

Z7T JO 

The form £) is a map from Lo 2n ®Lo 2n to the real numbers, and rj(6) and are 
elements of I/02«- The properties of the Killing form, as well as the decomposition 
of r]{9) and £(#) into Fourier modes, ensure that uo satisfies the conditions of a 
symplectic form, i.e. it is an antisymmetric, nondegenerate, closed two-form. Note 
that this equation defines a symplectic form for the loop group of any Lie group, 
with K standing for the Killing form of the corresponding Lie algebra. 

The Hamiltonian on S can be taken to be the same one as equation (fT5|) . and it 
is also a Hamiltonian on LC>2n via the embedding L0 2n C O res (H) (equation (|29|) ): 
we replace M(x,y) by M 7 (x, y) which is the integral kernel of M 1 (x) defined the 
same way as in (TT3"|) . To get equations of motion for mesons, we would normally 
need to derive the Poisson Brackets corresponding to the symplectic form given in 
equation fl35|) . However, we argue below that we can simply rely on the result for 
the U(N) case, i.e. that mesons on Gri satisfy 't Hooft's integral equation, and 
deduce that on the space S the equation for mesons will again be 't Hooft's integral 
equation, as is also true for QCD with an O(M) gauge group, to which we claim 
that the theory on S is equivalent: 

The relationship between L0 2n and S, i.e. the fact that L0 2n C O res (H) where 
S is the coset space O res {7i) /[0(TC+) x 0(TC-)) (equation ff25l) ). has an analog in 
the unitary case: it is true as well (see [7J) that LU n , the loop group of the unitary 
group, is a subgroup of U res (H), where we remember that Gri is the coset space 
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U res (?i)/[U(TC + ) x U(7i-)] (equation (TTU|) ). The symplectic form on Gri given 
in equation ffTTj) is invariant under the unitary group action on Gr% and is the 
unique such form (see pQ); similarly, the form in equation f|35|) . when it is taken 
on LU n , is invariant under the corresponding action in LU n (in fact, it is invariant 
under any translations in the group). Therefore, the meson equation derived on 
LU n using equation (1351) and the meson equation obtained on Gri using equation 
( TTTj) appear to be the same (or closely related), i.e they both are the 't Hooft 
integral equation. Similarly, equations obtained directly on S using a symplectic 
form invariant under the orthogonal group action, and the equations derived from 
equation fl35l) considered on LOm appear to be the same (or closely related). Since 
U res (Ti) and O res (TC) are analogous and the embeddings LU n C U res {7i) and L0 2n C 
O res (Ti.) have the same structure, the resulting meson equations on S and on LO<m 
are analogs of those on Gr\ and LU n , and we conclude that we get 't Hooft's 
equations on S (and on L02 n ) as we did on Gr\. 

Now that we established that our theory is consistent with 2DQCD through 
properties of mesons, it makes sense to consider the topological invariant of the 
space S, which is analogous to the virtual rank of Gri, to be related to baryon 
number, i.e. to consider baryons to be solitons, and we do so below. 

5 Baryon Number Non- Conservation 

As we have established, our phase space has two connected components (see equa- 
tions fl28l) and (l3Tj) ). and it is this topological property which we propose here to be 
related to the baryon number B. We will show in this section that the properties 
of baryon number given by our theory are consistent with QCD in two dimensions 
with an odd orthogonal gauge group. 
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We propose that in the space S, the topological invariant given in equation 
(|3ip corresponds to the parity of baryon number, i.e. states $ with = 

have even baryon number, and states with = 1 have odd baryon number. 

Correspondingly, in the space LOm^ shrinkable loops correspond to even baryon 
number, and non-shrinkable loops to odd baryon number. The physical meaning of 
relating a quantity to a topological invariant is that the quantity is conserved. So 
in this case, the quantity Qb defined by 

Q B = B mod 2 , (36) 

where B is the baryon number, is conserved, i.e. AQb = 0, which means that 
baryon number can change by multiples of two. In other words, baryon number is 
conserved only modulo two. Physically, this means that baryon-baryon pairs can 
annihilate 3 

The equivalence of the theory on S to QCD in two dimensions with an orthogonal 
gauge group is supported by the following argument: first, baryon-baryon annihila- 
tion can be shown directly to occur in O(M) QCD, where M here may be even or 
odd [9]: unlike SU(N), whose fundamental iV-dimensional representation is com- 
plex, O(M) has a real fundamental representation, which is the same as its dual. 
Since quarks are described by the fundamental representation while antiquarks are 
described by the dual of that representation, this implies that QCD with an O(M) 
gauge group does not distinguish quarks from antiquarks. Therefore, it does not 
distinguish baryons from antibaryons either, and just as a baryon-antibaryon pair 
can annihilate, so can a baryon-baryon pair. This is the same result we arrived at 
above! 

The idea that baryon number would be conserved only mod 2 appeared in a different context in [5]. 
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The question of baryon number conservation is a bit more subtle; baryon-baryon 
annihilation does not necessarily imply the conservation of baryon number modulo 
two: it means only that baryon number can be conserved at most modulo 2, i.e. 
single baryons might still appear or disappear, in which case there would be no 
conservation at all. We now show that baryon conservation is determined by the 
parity of M, i.e. it depends on whether the gauge group O(M) is an even orthogonal 
group, M = 2N, or an odd orthogonal group, M = 2N + 1. 

The lagrangian of QCD with an orthogonal gauge group has the same form 
as the one given in equation ([1]), except that the gauge fields are now the 
generators of the orthogonal rather than the unitary group (they are real anti- 
symmetric rather than complex anti-hermitian) and the quark field ip is real. It is 
clear that the lagrangian is invariant under the transformation ip — > —i/j, the only 
phase transformation allowed for real fields. This means that the only expectation 
values which can be non-zero are those which contain an even number of quarks, so 
that quarks are conserved modulo two. 

We can now see the difference between the even and odd cases. To form a baryon, 
we must take a product of M quarks. When M is even, the baryon is the product 
of an even number of quarks and therefore does not change under ip — > — ijj, which 
means that baryons are not conserved at all: Green's functions containing an odd 
number of baryons can be non-zero, and processes such as 2B — > 3B, which both 
begin and end with an even number of quarks, may occur. However, the situation 
when M is odd is different: each baryon consists of the product of an odd number 
of quarks so the baryon is odd under tp —>■ —ip. Therefore, non-zero correlation 
functions must necessarily contain an even number of baryons, which means that 
baryons are indeed conserved modulo 2, consistent with the theory on S. 
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6 Relation to master fields 



A classical field, such as the field M(x, y) which we have discussed, representing 
the large N limit of QCD is also known as a master field. There is another 
approach to master fields given in [JUJ [EE], where the mathematical formalism of 
non-commutative probability theory is used within a matrix model. A master field 
is in general a classical configuration such that in the large N limit, the values 
of gauge-invariant Green's functions are given simply by their value at the master 
field - i.e., no functional integral needs to be done. In [101 ITT] , the Hilbert space 
on which the master fields act is given by states generated from the vacuum by 
creation operators which satisfy the Cuntz algebra, i.e. 

a(x)\Vl >= 0, a(x)a'(y) = S(x — y), 

with no further relations. Master fields are given in terms of these creation and 
annihilation operators, i.e. M = M(a, a'). For the case of matrix fields which have 
independent distributions and are decoupled, there is the form 

M(a, a\ x) = a(x) + J2 M n+l a ]n (x) (37) 

n 

(the M n+ i turn out to be connected Green's functions). For the more general case 
in which the matrices are coupled, the form for the master field becomes 

oo 

M(x) = a(x) + ^*,yi,...,vfc at (2/i) ' ' ' a Kvk) ■ (38) 

k=l 

An interesting problem would be to make a connection between this framework 
and the master fields M(x, y) on the Grassmannian manifold described in Sections 
[3] and HJ i-e. to write M(x,y) in terms of these creation and annihilation operators. 
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A construction of the master field for 2DQCD is suggested in [10], but we propose 
a different form for M(a, a'): 

M{x,y) = a){x)a{y) . (39) 

This form is motivated by the intriguing fact that its commutation relations are 
identical to equation (I14p . except for the e term which does not have an analog here, 
and a sign: 

[M(x, y), M(z, w)} = 8(y - z)M(x, w) - 8(w - x)M(z, y) . (40) 

Therefore, it appears plausible that a matrix model with this master field could 
very well provide further insight into 2DQCD. A full investigation of this possibility 
is beyond the scope of this paper. 

7 Open Problems 

We have extended Rajeev's work to construct a classical dynamical system which is 
equivalent to the large N limit of two dimensional QCD with an 0(2N + 1) gauge 
group. We argued that the same equation of motion for mesons would result for the 
theory on S as for the theory on Gr±, which furthermore is the same equation derived 
within 2DQCD with 0(2N + 1) and U(N) gauge groups. Considering baryons as 
topological solitons in our theory, we showed that baryon number is conserved only 
modulo 2. We also showed that our model is related to the master field approach 
to matrix models by providing an explicit formula for the master field that has the 
same commutation relations as in our model. 

There are several directions for further investigations which may now be pursued. 
We suggest a few below: 
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It would be interesting to see what would happen when the symplectic groups 
or the exceptional groups serve as the QCD gauge groups. For Sp(2N), we can 
see immediately that a direct analog of our work would need to be modified. 
Suppose we propose an equivalence between a classical dynamical system on an 
infinite-dimensional Grassmannian manifold, or the loop group of the symplectic 
group, LSp2N, and the large N limit of 2DQCD with Sp(2N) gauge group. Since 
iri(Sp(2N)) = Z, our construction would lead to conservation of baryon number; 
however, baryons should not even exist in Sp(2N) QCD (they decay into mesons; 
see [9]). 

Another direction one might pursue is quantizing the classical dynamical system 
we have described, i.e. promoting the classical observables, which are just functions 
on the phase space, to operators. Commutation relations would be given by the 
Poisson Brackets. This has been done for the unitary case in [I], resulting in a 
theory of QCD at finite iV with a correspondence between h and 1/N. 

In addition, one may calculate the mass for a baryon by minimizing the Hamil- 
tonian subject to the constraint of odd baryon number. One would then interpret 
the resulting minimal energy configuration of odd baryon number to correspond to 
a stable state of one baryon. 

Another interesting investigation would be to study chiral symmetry breaking 
in this system and its relation to the non-conservation of baryon number. See [12] 
for a recent investigation of chiral symmetry breaking and stability of topological 
solitons in SU(N) and SO(N) Yang-Mills theories. 

Finally, a full-fledged investigation of the master field approach using the pro- 
posed master field (equation (139]) ) has the potential to shed new light on 2DQCD. 
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